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Abstract. A millimetre molecular line survey of seven high mass-loss rate carbon stars in both the northern and southern 
skies is presented. A total of 196 emission lines (47 transitions) from 24 molecular species were detected. The observed CO 
emission is used to determine mass-loss rates and the physical structure of the circumstellar envelope, such as the density and 
temperature structure, using a detailed radiative transfer analysis. This enables abundances for the remaining molecular species 
to be determined. The derived abundances generally vary between the sources by no more than a factor of five indicating that 
circumstellar envelopes around carbon stars with high mass-loss rates have similar chemical compositions. However, there are 
some notable exceptions. The most striking difference between the abundances are reflecting the spread in the '"C/'^C-ratio of 
about an order of magnitude between the sample stars, which mainly shows the results of nucleosynthesis. The abundance of 
SiO also shows a variation of more than an order of magnitude between the sources and is on the average more than an order of 
magnitude more abundant than predicted from photospheric chemistry in thermal equilibrium. The over-abundance of SiO is 
consistent with dynamical modelling of the stellar atmosphere and the inner parts of the wind where a pulsation-driven shock 
has passed. This scenario is possibly further substantiated by the relatively low amount of CS present in the envelopes. The 
chemistry occurring in the outer envelope is consistent with current photochemical models. 

Key words. Molecular processes - Stars: abundances - Stars: AGB and post-AGB - Stars: carbon - circumstellar matter 



1. Introduction 

The chemistry associated with carbon stars has long been 
known to be rich and complex in comparison to the alter- 
native O-rich regime (i.e., where C/0<1). This is in part 
due to the favourable bonding of the carbon atom, enabling 
long chains and complex species to form. Most of the cur- 
rent understanding of carbon stars has come from both ob- 
servational and theoretical work on the high-mass-losing car- 
bon star, IRC-i-10216. This source, which lies within 200pc 
and presents an ideal specimen for the study of carbon-rich en- 
velopes, has been mapped interferometrically in various molec- 



them. Much work has been done on the carbon-rich post-AGB 
sources CRL 618 and CRL 2688, and the chemistry of CRL 



ular species (e.g. 


Bieging & Tafalla 


1993; Dayal & Bieging 


1993 


, 1995|;|Gensheimeret al. 1995; 


guelin etal. 1993, 1996|; 


Lucas et al. 1995; 


Lucas & Guelin 


1999) and has had mod- 


els of its dust (a good summary is given by jMen'shchikov et al. 


2001 


) and chemistry (e.g. IVIillar et al 


200C|) constructed. These 



618 in particular has been modelled by Woods et al. (2002). 
Detailed chemical studies of carbon stars on the AGB have 
been few in number, but examples include the molecular line 
surve y of I RAS 15194-5115, a peculiar '■'C-rich star (Nyman 
et al. 1993 ). Carbon star surveys which include molecular-line 
comparisons ar e fewer, and have been limited in the number of 
lines observed. Olofsson et al. (1993a) detected some 40 stars 
in a h andful of species other than CO. The sample of Bujarrabal 
et al. ( 1994 ) included 16 carbon stars, with up t o ten molecular 
l ines o bserved in each. A more recent survey by [Olofsson et al 



(1998) detected 22 carbon stars in up to 6 molecular lines. 



tools have produced groundbreaking results and have been used 
to set a paradigm for what has come to be known as "carbon 
chemistry" in connection with evolved stars. However, the ac- 
curacy of employing IRC-i-10216 chemistry to similar carbon 
stars has been Uttle-tested due to the difficulties in observing 



The survey work presented here purports to be the most 
complete and consistent molecular-line survey in AGB carbon 
stars to date, covering high mass-loss rate objects in both the 
northern and the southern sky. Previously unpublished spectra 
of five stars (IRAS 15082^808, IRAS 07454-7112, CIT 6, 
AFGL 3068 and IRC-H40540) are presented, and spectra taken 
towards IRC+10216 and IRAS 15194-5115 with the Swedish- 
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ESQ Submillimetre Telescope (SEST; [Nyman et al.||1993D and 
IRC-H10216 with the Onsala Space Observatory (OSO) 20 m 
telescope are used to supplement the survey. Comparison of 
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Table 1. Positions, luminosities, periods and calculated distances of the sample of carbon stars. 



IRAS No. 


Other cat. name 


B1950 Co-ords. 


P 


L 


D 


T, 


Ti 


Li/L. 








[days] 


[Lo] 


[pc] 


[K] 


[K] 




07454-7112 


AFGL4078 


07:45:25.7 -71:12:18 




9 000° 


710 


1200 


710 


4.3 


09452 + 1330 


IRC+ 10216 


09:45:15.0 +13:30:45 


630 


9 600 


120 




510 




10131 +3049 


CIT6 


10:13:11.5 +30:49:17 


640 


9 700 


440 


1300 


510 


6.7 


15082 - 4808 


AFGL4211 


15:08:13.0 -48:08:43 




9 000" 


640 




590 




15194-5115 




15:19:26.9 -51:15:19 


580 


8 800 


600 


930 


480 


2.2 


23166+ 1655 


AFGL3068 


23.16.42.4 +16:55:10 


700 


7 800 


820 




1000 




23320 + 4316 


IRC+40540 


23:32:00.4 +43:16:17 


620 


9400 


630 


1 100 


610 


6.6 



(a) Assumed value. 



data from IRC+10216 taken with both telescopes affords a high 
degree of confidence in the relative calibration that can be de- 
rived. 

Up to 5 1 molecular lines were observed in the sample of 7 
high-mass-losing carbon stars, of which 47 were clearly de- 
tected. Mass-loss rates, dust properties and the '^CO/'^CO- 
ratio are calculated self-consistently using a radiative transfe r 
method ( [Schoier & Olofsson|^000|, |200l|;|Sch6ier et al.|^002| ). 
An approach similar to that of Nyman et al. ( 1993 ) is adopted 
to calculate fractional abundances (including upper limits), and 
a detailed analysis of the comparison between the calculated 
abundances is carried out. Hence, Sect. || details the observa- 
tions carried out and the instrumentation used. Section || gives 
the observational results, including a presentation of various 
spectra. Section ^ details the NLTE radiative transfer code used 
to determine the envelope parameters and Sect. |] explains the 
method of calculating chemical abundances. The results and 
deductions are discussed in Sect. ^. 

2. Observations 

2.1. Carbon star sample 



Following the CO survey of [Nyman et alj ( |1992| ), several car- 
bon stars, which were bright in CO lines, were selected for a 
more comprehensive molecular line search. These stars are rare 
in that they are all losing mass at a very high rate, and hence 
are more likely to produce strong emission from a variety of 
molecular lines. The sample of seven carbon stars (Table |l]) was 
observed u sing both the 15 m Swedish-ESO Submillimetre 
Telescope ( [Booth et al.l|l989| ) during the period 1987 - 1996, 
and the Onsala 20 m telescope, in 1994. The SEST, situated on 
La Silla, Chile, was used to observe IRAS 07454-7112, IRAS 
15082-4808, and IRAS 15194-5115. The 20m telescope, lo- 
cated at the Onsala Space Observatory (OSO) in Sweden, ob- 
served the remaining three sources, CIT 6, AFGL 3068 and 
IRC +40540. Both telescopes were used to observe the well- 
studied carbon star, IRC +10216 in order to determine the rel- 
ative calibration between the two telescopes. 

The JCMT public archive]] was searched for complemen- 
tary line observations, in particular the CO lines used in the 



radiation transfer modelling described in Sect. Lines for 
which multi-epoch observations are available in the JCMT 
archive typically display in tensities that are consistent to ~20% 
( Schoier & Olofsson 2001 ). In addition, interferometric obser- 
vations of the CO(7=1-0) brightness distr ibution around s ome 
of the sample stars have been performed (Nerietal 1998 ) us- 
ing the Plateu de Bure interferometer (PdBI), France. The data 
are publically available and have been used in this paper 

The coordinates used for each individual source are listed 
in Table |l} Also shown in Table |l] are the adopted luminosi- 
ties and distances to be used in the molecular excitation analy- 
sis. For stars where a period has be en determined (see Table p]) 
the period-luminosity relation from Groenewegen & Whitelock 
(1990 was used to estimate the corresponding luminosity. If a 
reliable period is not available the total bolometric luminosity 
was fixed to 9 000 L0. The distance was then obtained from the 
luminosity using the observed bolometric magnitude. Schoier 
& Olofsson (2001 ) used the same approach when determining 
the distances to a large sample of optically bright carbon stars 
and concluded that there were no apparent systematic effect 
when comparing with estimates based upon Hipparcos paral- 
laxes, although the scatter is large and the distance estimate for 
an individual source is subject to some uncertainty of up to a 
factor of ~2. The eff'ects of the adopted distance on the molec- 



ular excitation will be addressed in Sect. 5.5 



If, for simplicity, the central ratiation field is represented 
by one or two blackbodies their properties can be determined 
from a fit to the observed spectral energy distribution (SED) as 



described in Kerschbaum (1999 ). A fit to the SED gives the two 



blackbody temperatures (Tt and T^) as well as their relative 
luminosities (Ld/L»). The parameters obtained in this fashion 
are presented in Table |l[ 

2.2. Instrumentation 

The SEST is equipped with two acousto-optical spectrometers 
(HRS, 86 MHz bandwidth with 43 kHz channel separation and 
80 kHz resolution; LRS, 500 MHz bandwidth with 0.7 MHz 



' |ittp: //www. jch.hawaii.edu/3ACpublic/JCKT/ 
The JCMT is operated by the Joint Astronomy Centre in Hilo, 



Hawaii on behalf of the present organisations: the Particle Physics and 
Astronomy Research Council in the United Kingdom, the National 
Research Council of Canada and the Netherlands Organization for 
Scientific Research. 



Woods et al.: Molecular abundances in carbon-rich CSEs 



3 





SEST 15 m 


Onsala 20 m 


Frequency 


FWHM 77„,b 


FWHM 7,„,b 


[GHz] 


["] 


["] 



Table 2. Beam widths and efficiencies at selected frequencies. 1988). The '^CS (7=2-1) integrated intensities include both 

these lines. It is assumed that the relative intensities of these 
two lines are constant in all the objects, and hence will not af- 
fect comparisons. The lines of two HC3N isotopes, HC'^CCN 
and HCC'^CN, lie very closely and are in all cases blended 
together No attempt is made to separate the components. 

3.2. Line profiles 

There are four characteristic line profile shapes which give in- 
formation on the source being observed. These lineshapes are 
most typically seen in CO emission since it is often strongest 
(for a comprehensive review of line profiles see plofsson et al 



86 


57 


0.75 


44 


0.58 


100 


51 


0.73 


39 


0.53 


115 


45 


0.70 


33 


0.48 


230 


23 


0.50 






265 


21 


0.42 







channel separation and a resolution of 1.4 MHz). The receivers 
used were dual polarization Schottky receivers at both 3 and 
1.3 mm wavelength. Typical system temperatures above the at- 
mosphere were 400-500 K and 1000-1 800 K, respectively. 

The OSO 20 m telescope uses two filterbanks (MUL B, 
64 MHz bandwidth with a channel width of 250kHz; MUL A, 
5 12 MHz bandwidth and a channel width of 1 MHz). The re- 
ceiver used was a horizontally, linearly polarised SIS receiver 
with a typical system temperature of 400-500 K above the at- 
mosphere. 

All observations were performed using the dual beam 
switching method, which places the source alternately in two 
beams, and yields very flat baselines. Beam separation was 
in both cases about 11. '5. Calibration was done with the stan- 
dard chopper-wheel method. The intensity scales of the spec- 
tra are given in main-beam brightness temperature (the cor- 
rected antenna temperature, T^, divided by the main-beam effi- 
ciency, T/mb)- Main-beam efficiencies and FWHM beam widths 
are given in Table || for both telescopes. 

3. Observational results 

3.1. Observed lines 

A total of 196 lines were detected in the sample. 47 transitions 
of 24 molecular species were detected, and upper limits for an- 
other 95 relevant transitions were also obtained, including an- 
other three species (C2S, C3S and SO). Previously unpublished 
specti-a are shown in Figs. [l|-| & [a11 - |A.13[ Table | lists the 
detections in all seven sources, together with their peak and in- 
tegrated intensities. If a line has a hyperfine structure, the fre- 
quency and intensity of the strongest component is listed, and 
the integrated intensity is the sum over all hyperfine compo- 
nents. Values of T^b for lines where no detection was made are 
the rms noise values. 

Almost all lines observed in IRC -1-10216 by the SEST were 
observed in the same source with the OSO 20 m telescope. 
The majority of those which were not are due to the lack of a 
1.3 mm receiver at OSO. A large proportion of these lines were 
observed in the remaining six sources. The purpose of observ- 
ing IRC-i-10216 twice, with different telescopes, was to ascer- 
tain the relative calibration between the two different setups, 
and hence gain a basis from which good comparisons could be 
made between the entire sample of carbon stars. 

The '^CS (7=2-1) line observed in some of these sources 



1993b). For optically thick emission the line profile can be de- 
scribed as parabolic for unresolved sources or flat-topped for 
resolved sources. A parabolic profile is shown by most of the 
'^CO emission lines in the sample except in IRC H- 10216 and 
IRAS 15194-5115, where the CO (7=1-0) emission shows a 
flat-topped profile. 

When emission is optically thin, a flat-topped profile is seen 
for unresolved sources, and a double-peaked profile for a re- 
solved source. Examples of these profiles can be seen in the 
i^CO emission towards IRAS 15082-4808 and lRC+10216. In 
the following calculations, all lines are assumed optically thin, 
except those from CO and HCN. 

The values of the expansion velocity, Vgxp, of the circum- 
stellar envelopes quoted in Table ^ are obtained from the ra- 
diative transfer CO modelling of these sources (presented in 
Sect. ^ where its value is adjusted until a good fit to the ob- 
served line profiles is obtained. No trends in the widths of the 
observed lines are present which lends further support to the 
assumption adopted here that these envelopes are expanding at 
constant velocities. 



is partially blended with the C3S (7=16-15) line (ECahane et al 



Fig. 1. Low-resolution spectra of IRAS07454-71 12, obtained with the SEST. 



Fig. 2. Low-resolution spectra of IRAS07454-71 12, obtained with the SEST. 



Table 3. Detected lines 



Molecule 


Transition 


Frequency 


IRC+10216 (SEST) 


IRC+10216 (OSO) 


IRAS15194-5115 


IRAS 15082-4808 










J Tmbdv 


Tmb 


J Tmbdv 


Tmb 


J Tmbdv 


Tmb 


J Tnibdv 






[MHz] 


[K] 


[Kkms-1] 


[K] 


[Kkms-'] 


[K] 


[Kkms-'] 


[K] 


[Kkms-'] 


HC5N 


J=32-31 


85201.348 


0.15 


4.78 






<0.01 


<0.13 


<0.01 


<0.10 


C3H2 


2|,2-lo,i 


85338.905 


0.11 


2.49 






0.02 


0.84 


<0.01 


<0.10 


C4H 


N=9-8, J=19/2-17/2 


85634.00 


0.20 


9.60 










<0.01 


< 0.05 


C4H 


N=9-8, J=17/2-15/2 


85672.57 


0.20 


9.60 










<0.01 


< 0.05 


C2S 


6(7)-5(6) 


86181.413 


<0.03 


<0.46 


<0.04 


<0.76 


<0.07 


<0.11 


< 0.04 


< 0.06 


h"cn 


J=l-0 


86340.184 


3.23 


87.20 


3.86 


114.14 


0.57 


20.93 


0.09 


3.06 


C3S 


J=15-14 


86708.379 


< 0.02 


< 0.30 


< 0.03 


<0.65 


< 0.01 


< 0.08 


<0.01 


< 0.06 


SiO 


J=2-l 


86846.998 


0.71 


20.40 


0.95 


24.48 


0.16 


5.33 


0.07 


1.84 


hn"c 


J=l-0 


87090.859 


< 0.04 


< 0.62 






0.02 


1.11 


<0.01 


<0.11 


C2H 


N=l-0 


87316.925 


0.52 


41.60 


0.66 


49.48 


0.13 


11.50 


0.04 


4.04 


HCN 


J=l-0 


88631.847 


7.80 


190.93 


9.14 


214.04 


0.55 


15.33 


0.37 


11.26 


C3N 


N=9-8, J=19/2-17/2 


89045.59 


0.31 


9.19 








<0.01 


<0.07 


<0.01 


<0.09 


C3N 


N=9-8, J=17/2-15/2 


89064.36 


0.31 


9.19 








<0.01 


<0.07 


<0.01 


<0.09 


HC5N 


J=34-33 


90525.892 


0.16 


2.58 


<0.03 


<0.67 








<0.01 


<0.10 


HC"CCN 


J=10-9 


90593.059 


0.09 


1.98 


<0.03 


<0.34 


0.02 


0.76 


<0.01 


<0.05 


HCC^CN 


J=10-9 


90601.791 


0.09 


1.98 


<0.03 


<0.34 


0.02 


0.76 


<0.01 


<0.05 


HNC 


J=l-0 


90663.543 


0.74 


26.49 


0.58 


17.22 


0.08 


3.72 


0.05 


1.70 


SiS 


J=5^ 


90771.546 


1.14 


35.54 


1.07 


39.47 


0.05 


2.78 


0.03 


1.58 


HC3N 


J=10-9 


90978.993 


2.42 


64.59 


2.00 


52.63 


0.09 


3.71 


0.11 


4.15 


"cs 


J=2-l 


92494.299 


0.07 


2.97 








0.06 


2.67 


<0.01 


<0.10 


SiC2 


4o,4-3o,3 


93063.639 


0.46 


14.60 














0.03 


0.40 


HC5N 


J=35-34 


93188.127 


0.12 


5.50 




^ 


_ 





<0.01 


<0.10 


C4H 


N=10-9, J=21/2-19/2 


95150.32 


0.19 


6.08 


0.18 


8.75 


0.04 


1.95 


0.02 


0.47 


C4H 


N=10-9, J=19/2-17/2 


95188.94 


0.19 


6.08 


0.18 


8.75 


0.04 


1.95 


0.02 


0.47 


SiC2 


42,2-32,1 


95579.389 


0.32 


10.73 


— 


— 


— 


— 


<0.01 


<0.08 


C^*S 


J=2-l 


96412.982 


0.15 


5.66 


— 


— 


— 


— 


— 


— 


cs 


J=2-l 


97980.968 


2.92 


81.23 


4.06 


118.89 


0.37 


15.21 


0.25 


8.44 


C3H 


2ni/2,9/2-7/2 


97995.450 


0.18 


11.27 


0.28 


14.54 


0.01 


0.64 


<0.01 


<0.11 


so 


J=3-2 


99299.879 


<0.02 


<0.25 




















HC3N 


J= 11-10 


100076.389 


2.21 


53.84 


— 


— 


— 


— 


— 


— 


HC"CCN 


J=12-ll 


108710.523 


0.09 


2.34 


0.14 


1.22 


0.01 


0.39 


<0.01 


<0.03 


HCC"CN 


J=12-ll 


108721.008 


0.09 


2.34 


0.14 


1.22 


0.01 


0.39 


<0.01 


<0.03 


"CN 


N=l-0 


108780.201 


0.11 


8.56 


0.24 


6.07 


0.02 


2.18 


<0.01 


<0.05 


C3N 


N=ll-10, J=23/2-21/2 


108834.27 


0.39 


11.76 


0.64 


18.00 


0.01 


0.68 


0.02 


0.63 


C3N 


N=ll-10, J=21/2-19/2 


108853.02 


0.39 


11.76 


0.64 


18.00 


0.01 


0.68 


0.02 


0.63 


SiS 


J=6-5 


108924.267 


1.10 


30.00 


1.86 


52.40 


0.05 


2.44 


0.04 


1.14 


HC3N 


J=12-ll 


109173.634 


2.03 


52.39 








0.07 


3.16 


0.08 


2.79 


c«o 


J=l-0 


109782.160 


0.02 


0.76 














<0.01 


<0.17 


"CO 


J=l-0 


110201.353 


1.69 


26.62 


2.50 


36.40 


0.31 


14.93 


0.06 


1.82 


CH3CN 


6(1)-5(1) 


110381.404 


0.09 


4.77 


<0.05 


<1.42 


<0.01 


<0.18 


<0.01 


<0.23 


C2S 


8(9)-7(8) 


113410.207 


<0.03 


<0.48 


<0.07 


<1.96 


<0.02 


<0.53 


<0.01 


<0.11 


CN 


N=l-0 


113490.982 


3.43 


238.43 


3.88 


310.20 


0.13 


16.43 


0.19 


14.71 


CO 


J=l-0 


115271.204 


10.29 


269.43 


20.83 


542.71 


1.27 


55.86 


1.11 


36.71 


SiCj 


5o,5-4o,4 


115382.38 


0.80 


21.43 


1.60 


39.38 


0.07 


2.66 


0.03 


0.82 


c'^o 


J=2-l 


219560.319 


0.15 


4.21 














"CO 


J=2-l 


220398.686 


3.85 


79.62 






0.90 


35.58 


0.15 


5.76 


CH3CN 


12(0)-11(0) 


220747.268 


0.17 


2.80 






<0.05 


<1.01 


<0.03 


<0.48 


CN 


N=2-l 


226874.564 


2.55 


231.18 






0.15 


17.06 


0.20 


11.32 


CO 


J=2-l 


230538.000 


34.60 


799.00 






4.20 


150.00 


2.54 


69.40 


CS 


J=5^ 


244935.606 


16.57 


405.53 






0.77 


32.34 


0.47 


10.89 


HCN 


J=3-2 


265886.432 


45.10 


1010.24 






4.69 


139.76 


0.88 


20.83 



Table 3. (cont.) Detected lines 



Molecule 


Transition 


Frequency 


IRAS07454-7112 




C1T6 


AFGL 3068 


IRC+40540 








7"mb 


J Tm\jdv 


Tmb 


J Tmbdv 




J Tmbdv 




J T^mbdv 






(MHz) 


(K) 


(Kkms"') 


(K) 


(Kkms-1) 


(K) 


(Kkms"') 


(K) 


(Kkrns"') 


C2S 


6(7)-5(6) 


86181.413 


< 0.01 


<0.07 


<0.02 


<0.36 


<0.02 


<0.44 


<0.01 


< 0.28 


h"cn 


J=l-0 


86340.184 


07 


1.45 


0.15 


3.13 


0.13 


3.76 


0.13 


3.85 


C3S 


J=15-14 


86708.379 


<0.01 


<0.07 


<0.01 


<0.23 


<0.02 


< 0.47 


<0.01 


<0.21 


SiO 


J=2-l 


86846.998 


0.03 


0.52 


0.17 


3.16 


<0.02 


<0.47 


0.04 


0.94 


hn"c 


J=l-0 


87090.859 


< 0.01 


<0.08 


<0.01 


<0.23 


< 0.02 


< 0.47 


< 0.01 


< 0.25 


C2H 


N=l-0 


87316.925 


< 0.01 


< 0.09 


0.10 


5.81 


0.09 


5.94 


< 0.01 


< 1.78 


HCN 


J=l-0 


88631.847 


0.19 


4.43 


0.77 


18.07 


0.39 


8.97 


0.26 


8.24 


C3N 


N=9-8, J= 19/2-17/2 


89045.59 


<0.01 


<0.07 


— 


— 










C3N 


N=9-8, J= 17/2-15/2 


89064.36 


<0.01 


<0.06 


— 


— 










HC5N 


J=34-33 


90525.892 


0.02 


0.06 


<0.01 


<0.21 


<0.01 


<0.21 


<0.02 


<0.37 


HC"CCN 


J=10-9 


90593.059 


0.01 


0.13 


<0.01 


<0.11 


<0.01 


<0.10 


0.08 


0.35 


HCC"CN 


J=10-9 


90601.791 


0.01 


0.13 


<0.01 


<0.11 


<0.01 


< 0.10 


0.08 


0.35 


HNC 


J=l-0 


90663.543 


0.02 


0.50 


0.13 


2.95 


0.06 


1.18 


0.09 


1.42 


SiS 


J=5^ 


90771.546 


0.03 


1.07 


0.04 


0.59 


0.04 


0.59 


0.15 


4.11 


HC3N 


J=10-9 


90978.993 


0.06 


1.39 


0.19 


5.17 


0.12 


3.37 


0.19 


6.16 


"cs 


J=2-l 


92494.299 


<0.01 


<0.04 
















SiC2 


4o,4-3o,3 


93063.639 


0.01 


0.13 


— 


— 










HC5N 


J=35-34 


93188.127 


<0.01 


<0.04 


— 


— 










C4H 


N=10-9, J=21/2-19/2 


95150.32 


<0.01 


<0.03 


<0.01 


<0.13 










C4H 


N=10-9, J=19/2-17/2 


95188.94 


<0.01 


<0.03 


<0.01 


<0.13 










SiC2 


42,2-32,1 


95579.389 


<0.01 


<0.06 


— 


— 










cs 


J=2-l 


97980.968 


0.13 


2.53 


0.67 


16.91 


0.14 


3.93 


0.30 


9.56 


C3H 


2ni/2,9/2-7/2 


97995.450 


<0.01 


<0.09 


<0.02 


<0.55 


<0.01 


<0.31 


<0.02 


<0.45 


HC"CCN 


J=12-ll 


108710.523 


0.01 


0.08 


0.05 


0.44 


<0.01 


<0.14 


<0.01 


<0.09 


HCC"CN 


J=12-ll 


108721.008 


0.01 


0.08 


0.05 


0.44 


<0.01 


<0.14 


<0.01 


<0.09 


"CN 


N=l-0 


108780.201 


yJ.yJL 


ft 7'^ 


0.05 


4.21 


^ ft fti 


<. u.zo 


^ n ni 

<^ U.Ul 


^ ft 1 R 


C3N 


N= 11-10, J=23/2-21/2 


108834.27 




ft O'i 


0.10 


2.19 


ft ft*\ 




U.Uj 


ft 70 

u. /u 


C3N 


N=ll-10, J=21/2-19/2 


108853.02 


0.01 


0.23 


0.10 


2.19 


0.05 


1.34 


0.03 


0.70 


SiS 


J=6-5 


108924.267 


0.03 


0.88 


0.19 


3.62 


0.15 


2.23 


0.07 


1.91 


HC3N 


J=12-ll 


109173.634 


0.06 


1.48 
















"CO 


J=l-0 


110201.353 


ft ftO 


2 10 


0.40 


5.90 


ft 17 


y.jj 


22 




CH3CN 


6(1)-5(1) 


110381.404 


^'ft ftl 


<0 12 


<0.01 


<0.29 


^'ft ft'^ 


^ U.OJ 


<. U.Uj 


ft 79 


C2S 


8(9)~7(8) 


113410.207 


<0.01 


<0.10 


<0.03 


<0.89 


<0.02 


< 0.53 


< 0.02 


< 0.64 


CN 


N=l-0 


113490.982 


0.17 


10.30 


1.12 


61.84 


0.09 


5.97 


0.41 


21.84 


CO 


J=l-0 


115271.204 


1.16 


25.57 


3.52 


105.00 


2.08 


47.08 


1.94 


41.88 


SiC2 


5o,5^0,4 


115382.38 


<0.01 


<0.24 


0.23 


6.71 


<0.02 


<0.62 


<0.03 


<0.75 


"CO 


J=2-l 


220398.686 


0.25 


7.34 














CH3CN 


12(0)-11(0) 


220747.268 


0.04 


0.67 














CN 


N=2-l 


226874.564 


0.41 


20.59 














CO 


J=2-l 


230538.000 


2.60 


50.00 














CS 


J=5^ 


244935.606 


0.49 


8.12 














HCN 


J=3-2 


265886.432 


0.50 


8.82 
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3.3. Upper limits 

The integrated intensities for lines which were not detected are 
determined using. 



•'exp 



Av = 3o-{2ny'^Av. 



(1) 



where cr is the rms noise in the spectrum, Av is the frequency 
resolution of the spectrum, Vjes the velocity resolution of the 
spectrum and n the number of channels covering the line width. 
Integrated intensities calculated with this method were used to 
determine upper limits to abundances. 

4. Full radiative transfer modelling 

In order to obtain accurate values of the mass-loss rates for 
the objects in the sample, detailed non-LTE radiative trans- 
fer modelling of the CO line emission was carried out. This 
modelling has proven to be one of the most reliable meth- 
ods for estimating many of the parameters characterising a 
circumstellar envelope, in particular, the mass-loss rate, ki- 



netic temperature structure, and expansion veloci ty ( [Kastner 
' 1992| ; |Groenewegen||l994|;|Crosas & Menterj|l997|; Sch oier & 



Olofsson |200l| ; |Schoieretal.||2002| ; |Olofsson et al.||2002| ). Once 
the basic envelope parameters are known (as well as the stel- 
lar parameters), the abundances of various molecules present 
in the circumstellar envelopes can be determined. 

4.1. CO line modelling 



Here the detailed non-LTE radiative transfer code of [Schoier 



(2000) is adopted and the modelling procedure as outlined in 



detail in Schoier & Olofsson (2001 ) is used. The code has been 



tested against a wide variety of molecular line radiative trans- 
fer codes, for a number of benchmark problems, to a high ac- 



( Bernes 



curacy (van Zadelhoff et al.| |2002[ ). The Monte Carlo method 



1979) is used to determine the steady-state level popu- 



lations of the CO molecules in the envelope as a function of dis- 
tance from the star, using the statistical equilibrium equations. 
In addition, the code simultaneously solves the energy balance 
equation including the most relevant heating and cooling pro- 
cesses. Heating is dominated by collisions between the dust 
and gas except in the outermost parts of the envelope where 
the photoelectric effect effectively heats the envelope. Cooling 
is generally dominated by molecular line cooling from CO but 
adiabatic cooling due to the expansion of the envelope is also 
important. The excitation analysis allows for a self-consistent 
treatment of CO line cooling. 

The observed circumstellar CO line emission is modelled 
taking into account 50 rotational levels in each of the funda- 
mental and first excited vibrational states. The energy levels 
and radiative transition probabilities from Chandra et al. (1996) 
are used. The recently published coUisional rates of CO with 
H2 by Flowei (2001) have been adopted assuming an ortho-to- 



para ratio of 3. For temperatures above 400 K the rates from 



Schinke et al, (1985) were used and further extrapolated to 



include transitions up to J-50. The collisional rates adopted 
here differ from those used in the previous modelling of some 



of these sources dRyde et al.|[l999| [Schoier & Olofsson) |2000i 
2001) and account for the slightly different envelope parame- 
ters derived in the present analysis. The same set of collisional 
rates were used for all CO isotopomers. 

The envelopes are assumed to be spherically symmetric 
and to expand at a constant velocity and the model includes 
the radiation emitted from the central star Dust present around 
the star will absorb parts of the stellar radiation and re-emit 
it at longer wavelenghts. For simplicity, the central radiation 
field is represented by one or two blackbodies and is deter- 
mined from a fit to the observed spectral energy distribution 
(SED) (Table 1). The inner radius of the circumstellar envelope 
is taken to reside outside that of the central blackbodies. This 
procedure provides a good description of the radiation field to 
which the envelope is subjected. For the sample stars, which 
have dense CSEs, the line intensities derived from the '^CO 
model are not sensitive to the adopted description of the stellar 
spectrum due to the high line optical depths. The stellar photons 
are typically absorbed within the first few shells in the model. 
A more detailed treatment of thermal emission from the dust 
present in the CSE, and the increase of total optical depth at 
the line wavelenghts, has been found to be of no major impor- 
tance in deriving the envelope parameters for high mass-loss 
rate objects (Schoier et al. 2002). 

The abundance of '^ CO relative to H2 was fixed at 1.0 x 
10"^, in a greement with Willacy & Cherchneff ( 1998 ) and the 
survey of Olofsson et al. ( 1993a ). The CO env elope size was 
estimated based upon modelling results from [Mamon et al 



(1988), which have been shown to compare well with obser- 
vations ( [Schoier & Olofsson] ^00l| ). The same envelope size 
was assumed for all CO isotopomers. 



4.2. Fitting strategy 

With the assumptions made in the standard model there re- 
mains two principal free parameters, the mass-loss rate (M) and 
the so-called /z-parameter, which contains information about in- 
dividual dust grains and controls the amount of heating through 
dust-gas collisions. 



p.O g cm-3 \ / 0.05 yum 

ooT 



Pd 



ad 



(2) 



where T is the dust-to-gas mass-loss rate ratio, pa the mass 
density of a dust grain, and oa the size of a dust grain, and the 
adopted values of these parameters are also given. The mass- 
loss rate and the /i-parameter were allowed to vary simultane- 
ously. Once the molecular excitation, i.e., the level populations, 
is obtained the radiative transfer equation can be solved exactly. 
The resulting brightness distribution is then convolved with the 
appropriate beam to allow a direct comparison with observa- 
tions. The beam profile used in the convolution of the modelled 
emission is assumed to be Gaussian which is appropriate at the 
frequencies used here. The best fit circumstellar model for a 
particular source is estimated from a;\f^-analysis (for details see 
[Schoier & Olofsson] pOOTj ) using observed '^CO integrated in- 
tensities and assuming a 15-20% calibration uncertainty. The 
number of observational constraints used in the modelling is 
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Table 4. Summary of circumstellar properties derived from the CO modelling (see text for details) 
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Source 


M 


h 




^red 




(M/Vcxp)* 


4 
p 


d 


'^CC 




[Moyr-'] 








[km s"'] 


( M /Wvn)Tn91f^ 
\ / CXU / 1 UZ ID 


[cm] 


[cm] 


"CO 


IRAS07454-7112 


5.0 (-6) 


1.5" 


2 


0.1 


13.0 


0.46 


3.0(17) 


2.4(16) 


17 


IRC+10216 


1.2 (-5) 


1.5 


27 


1.4 


14.5 


1.00 


3.7(17) 


5.5(16) 


45 


CIT6 


5.0 (-6) 


2.3 


20 


0.8 


17.0 


0.36 


1.9(17) 


2.8(16) 


35 


IRAS 15082-4808 


1.0 (-5) 


1.5" 


2 


0.2 


19.5 


0.62 


2.5(17) 


3.5(16) 


35 


IRAS 15 194-5 115 


1.2 (-5) 


2.8 


10 


0.9 


21.5 


0.67 


3.2(17) 


7.2(16) 


6 


AFGL3068 


2.0 (-5) 


2.5 


4 


1.6 


14.0 


1.72 


3.8(17) 


1.2(17) 


30 


IRC +40540 


1.5 (-5) 


1.5 


7 


1.0 


14.0 


1.29 


4.0(17) 


7.3(16) 


50 



(a) Adopted value 

(b) Number of observational constraints used in the '^CO modeUing 

(c) Reduced of the best fit '^CO model 

(d) Tp is the photodissociation radius of CO 

(e) Determined from radiative transfer modelling of both '"CO and ''CO emission 



shown in Table ^ together with the reduced of the best fit 
model. In all the cases indicating good fits. 

The '^CO data used in the analysis are presented in Schoier 



(2002) for CIT 6, 



& Olofsson ( pOOlD and [Schoier et al 
IRC+10216, IRAS 15194-5115, and IRC+40540 and consist 
of both millimetre and sub-millimetre line data as well as far 
infra-red high-7 transitions observed by ISO. The best fit '^CO 
model obtained for AFGL 3068 is overlayed onto observa- 
tions and presented in Fig. |. For IRAS 07454-71 12 and IRAS 
15082^808 only 7=1-0 and 7=2-1 Une data as observed by 
the SEST, and presented here in Figs. [l|, |, ^ and |a!6| - 
|A.9| , are available. Due to the limited number of constraints 
the /i-parameter was assumed to be equal to 1.5, i.e., the dust 
properties were taken to be the same as for IRC+10216 and 
IRC+40540 for these two sources. 



4.3. The mass-loss rates 

The determination of accurate mass-loss rates for some of 
these carbon stars is, it turns out, difficult in models where 



CO cooling is treated in a self-consistent manner ( [Sahai| |1990| ; 
Kastnei 1992; Schoier & Olofsson 2001). The reason for this 



is that an increase in M leads to more net cooling than heat- 
ing, which compensates for the increase in molecular density, 
making it hard to simultaneously constrain both M and the 
/i-parameter For the sources CIT 6, IRC+10216 and IRAS 
15194-51 15, where a significant number of high-7 transitions 
were observed by ISO, the degeneracy between M and the h- 
parameter is partly lifted allowing for better constraints to be 
put on these parameters (Schoier et al. 2002). For these three 
sources the mass-loss rate is estimated to be accurate to about 
50% within the adopted circumstellar model. For the remain- 
ing sources the uncertainty in the derived mass-loss rate is 
larger due to either the above mentioned degeneracy (AFGL 
3068 and IRC+40540) or that only two '^CO lines are used 
in the analysis (IRAS 07454-7112 and IRAS 15082-4808). 
Radiative transfer analysis of the observed continuum emis- 
sion from these sources (Schoier et al., in prep) give mass-loss 
rates that agree within ~50% when compared with those de- 



rived from the CO modelling. This is reassuring and lends fur- 
ther credibility to the mass-loss rates presented in Table ^. 

The CSEs of the sample stars have similar physical prop- 
erties. However, the stars presented here are losing mass at a 
significantly higher rate than the average carbon star: Schoier 
& Olofsson ( 2001 ) measure a median mass-loss rate for carbon 
stars of ~3xlO"^ Mq yr"', based on a sample of carbon stars 
complete within ~ 600 pc from the sun. This suggests that the 
stars in the sample presented here are going through the super- 
wind phase of evolution (e.g. Vassiliadis & Wood 1993 ) at the 
end of the AGB, and will soon eject the entire stellar mantle. 

4.4. CO isotopic ratios 

Using the envelope parameters derived from the '^CO mod- 
elling, the abundance of '^CO and C'**0 was estimated using 
the same radiative transfer code. The derived '^CO/'^CO-ratios 
are presented in Table Q The best fit '^CO model obtained for 
AFGL 3068 is overlayed onto observations and presented in 
Fig. H The fits to the observed '^'CO emission are good in all 
cases withx^^^^^- The quality of the fits further strengthen the 
adopted physical models for the envelopes. The '^CO/'^CO- 
ratio derived in this manner allows the assumption of optically 
thin emission adopted in Sect. ^ to be tested for those molecu- 
lar species where the isotopomer containing '^C has been de- 
tected. All the sample stars, except for IRAS 15194-5115 (to 
be discussed in Sect. 5.2), have inferred '^C/'^C-ratios in the 
range ~ 20 - 50. For optically bright carbon stars, i.e. gener- 
ally lower mass-loss rate objects, values between ~ 20-90 



are found (e.g., [Schoier & Olofsson 2000; Abia et al. 2001). 
The observed '^C/'^*C-ratios are in agreement with evolution- 
ary scenarios where the '^C/'^C-ratio is thought to increase 
from an initial low value of ~ 10 - 25, that depends on stellar 



mass, as the star evolves along the AGB (Abia et al. 2001 ). The 



thermal pulses that an AGB star experiences will effectively 
dredge nucleosynthesized '^C to the surface (Iben & Renzini 



1983) and, in addition to increasing the '^C/'-'C-ratio, eventu- 



ally turn it into a carbon star with C/0>1 in its atmosphere. 

C'^O emission was only detected towards IRC+10216 and 
a C'^O/C'^O-ratio of 1050 is derived, using four observational 
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-50 _ -50 _ -50 _ 10 20 30 

V [km s ^] V [km s ^] v [km s ^] Offset ["] 



Fig. 3. Multi-transition CO millimetre-wave line emission observed towards AFGL 3068. The observed spectra (histograms) 
have been overlayed with the best fit model results (solid lines). Also shown, lower right panel, is the observed radial brightness 
distribution (boxes with error bars) overlayed by the results from the model (full line), with the circular beam used in the 
radiative transfer analysis (dot-dashed line). The transition, telescope used, and the corresponding beamsize, are indicated for 
each observation. 



constraints (including three different transitions J -1-0, 2-1, 
3-2) for C'**0. This value is in e xcellent agreement w ith the 
i60/'^0-ratio of 1260 obtained by |Kahane etall ( |l992| ), using 
a combination of optically thin emission lines. In comparison, 
the value of this ratio in the solar neighbourhood is around 500. 
The '^O/'^O-ratio and in particular the '^0/"^0-ratio, which is 
not measured here, can be used as tracers of nucleosynthesis. 
Like the '^C/'-'C-ratio, these rati os are thought to inc rease as 
the star evolves along the AGB. Kahane et al. (1992) indeed 
find support for this scenario in a small sample of carbon-rich 
evolved stars. 



5. Abundance calculations 

A full radiative transfer analysis of the wealth of molecular data 
presented here is beyond the scope of this paper. A detailed 
treatment of the excitation would have to include the effects 
of dust emission and absorption in addition to accurate rates 
for collisional excitation of the molecules. Instead, abundances 
have been calculated for species with emission that is expected 
to be optically thin. In the case that the emission is optically 
thick, the calculated abundances using this assumption will be 
lower limits. 

The calculation of isotope ratios of various species (as 
shown in Table H) shows in general that many lines are optically 
thick (i.e., the ratio derived from observations is lower than the 
'^CO/'^CO abundance ratio derived from the radiative transfer 
analysis). Where this is the case, the abundance of the main iso- 
tope (viz., CN, CS and HC3N) has been calculated by scaling 
the abundance of the less abundant isotope by the '^CO/'^CO- 
ratio. This is indicated by the bold-faced type in Table ^ The 



abundance of HCN is not calculated using Eq. ^ since the line 
is certainly optically thick, but in all cases only by scaling the 
abundance of H'^CN by the calculated '^CO/'^CO-ratio. 

5.1. Radiative transfer 

In the calculation of molecular abundances the circumstellar 
envelope is assumed to have been formed by a constant mass- 
loss rate and to expand with a constant velocity, such that the 
total density distribution follows an r^^ law. It is further as- 
sumed that the fractional abundance of a species is constant in 
the radial range r, to and zero outside it. The excitation tem- 
perature was assumed to be co nstant throughout the emitting 
region. With these assumptions Olofsson et al. ( 199C ) showed 
that for a given molecular transition. 



Tmbdv 



2kv^ 

guAul (1 



[Bv(rex) - Bv(rbg)] X 



-41n(2)j:^ 



Ax, 



(3) 



where c is the speed of light, h the Planck constant, k the 
Boltzmann constant. By the Planck function, T^^ the excitation 
temperature, Tbg is taken to be the blackbody temperature of 
the cosmic background radiation at 2.7 K, mn, is the mass of 
an H2 molecule, Vgxp the expansion velocity of the circumstellar 
envelope, B the FWHM of the telescope beam, D the assumed 
distance, Mh, the mass-loss rate, Z the partition function, v the 
frequency of the line, gu the statistical weight of the upper level, 
Aui the Einstein coefficient for the transition, E\ the energy of 
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Table 6. Photodissociation rates 



r [cm] 



Fig. 4. The density (solid line) and kinetic temperature (dotted 
line) structures obtained from the CO excitation analysis for 
AFGL 3068. Also shown is the excitation temperature of the 
CO(7=2-l) hne (dash-dotted line). 



the lower transition level, and jci e - ri e/BD. The integral over 
X takes care of the beam filling. However, if r <s; BD, the in- 
tegral can be simplified to {r^ - r\)IBD. This is the case in the 
furthest sources, IRAS 07454-7112, CIT 6 and IRAS 15082- 
4808. However, for the CN emission from these sources, which 
tends to be very extended, and for the remaining four sources, 
the full integral is calculated. 



5.2. The excitation temperature 

Of importance in the radiative transfer analysis is the ex- 
citation (rotational) temperature assumed for the molecular 
emission. In Fig. ^ the kinetic temperature and density struc- 
tures for AFGL 3068 are shown. In the region where most 
of the molecular emission observed is thought to emanate 
from, ~ lO'^ - 10'^ cm, the kinetic temperature ranges from 
~ 30 - 10 K. Thus, if the emission were to be fully thermalized, 
then excitation temperatures would be expected to lie within 
this range. However, in addition to the low temperatures, the 
relatively low densities in this region, ~ 10^ - 10^ cm"^, makes 
the CO emission partly sub-thermally excited as illustrated in 
Fig. ^ Thus, sub-thermal emission is to be expected for all 
molecular species and the excitation temperature will depend 
on the species (and transition) observed. 

When two or more transitions of the same molecule are ob- 
served, it is possible to make an estimate of the rotation temper- 
ature (Tiot) of that molecule using Eq. ^ Assuming a molecular 
species to be excited over the same radial range, and according 
to a single temperature, the rotational temperature can be esti- 
mated. The results are shown in Table ^ It is clear that the ro- 
tational temperatures vary from source to source and between 
molecular species in the range ~ 3 - 30 K, as to be expected. 
The average excitation temperature is 8.7 K (averaged over the 



Molecule 


Go 
[s-'] 




HCN 


1.1 10^ 


1) 


CN 


2.5 10- 


10 


C2H2 


2.1 10- 


9 


CjH 


3.4 10- 


10 


SiS 


6.3 10- 


10 


cs 


6.3 10- 


10 


SiO 


6.3 10- 


10 



individual excitation temperatures for all sources, rather than 
molecular species). A generic value of 10 K was assumed for 
all the abundances estimated. 



5.3. The partition function 

All molecules are assumed to be linear, rigid rotators, except 
for SiCi and C3H2 which are asymmetric tops and CH3CN, 
which is a prolate symmetric top. Einstein A-coefficients 
(where a vailab le) and energy le vels are taken from Chandra & 
Rashmi ( |l998|) for SiCa from [Vrtilek et al.| ( |l987| ) for C3H2, 
and from [Boucher et al.| ( |l98C| ) for CH3CN. 

The partition function, Z, is calculated assuming that no 
molecules have a hyperfine structure, i.e. as having simple rota- 
tional energy diagrams. Hence the integrated intensities which 
are summed over all hyperfine components are used. For C3H2, 
CH3CN and SiC2 we use the approximate exp ression for an 
asymmetric rotor ( Townes & Schawlow 1975 ) multiplied by 
2 for C3H2, by 4 for CH3CN (to account for spin statistics) 
and by 1/2 for SiC2 (since half of the energy levels are miss- 
ing because of spin statistics). It is assumed that all levels are 
populated according to Tgx. 



5.4. Sizes of emission regions 
5.4.1. Photodissociation radii 

The chemical richness observed towards carbon stars can 
be qualitatively understood in terms of a photodissociation 
model (see the review by Glassgold 1996). Carbon-bearing 
molecules like CO, C2H2, HCN and CS, in addition to Si- 
bearing molecules like SiS and SiO, are all thought to be of 
photospheric origin where they are formed in conditions close 
to LTE. The photodissociation of these so-called parent species 
produces various radicals and ions that, in turn, drive a com- 
plex chemistry through ion-molecule and neutral-neutral reac- 
tions. For example, the radicals C2H and CN are the photodis- 
sociation products of C2H2 and HCN, respectively. In addition, 
the formation of dust in the CSE will affect the abundances of 
some of the species, in particular SiO and SiS, which condense 
onto dust grains. All other species observed in the sample are 
thought to be products of the circumstellar chemistry. 

To calculate the radial extent of the molecules HCN, 
CN, C2H and CS the photodissociation model of Huggins & 
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Table 5. Rotation temperatures 



Source 


"CO 


CN 


CS 


C3N 


C4H 


HCN 




2-1/1-0 


2-1/1-0 


5^/2-1 


11-10/10-9 


10-9/9-8 


3-2/2-1 


IRAS 07454-7112 


7.0 


5.2 


11.3 






5.7 


IRC+ 10216 (S) 


6.3 


3.9 


13.7 


17.9 


5.3 


7.8 


IRAS 15082-4808 


6.5 


3.6 


8.3 






5.6 


IRAS15194-5115 


5.6 


4.0 


9.7 






9.7 


Average 


6.4 


4.2 


10.8 


17.9 


5.3 


7.2 




Source 


H"3'C3N 


HC3N 


HC5N 


SiC2 


SiS 






12-11/10-9 


12-11/10-9 


34-33/32-31 


5-4/4-3 


6-5/5^ 


Average 


1RAS07454-7112 


8.5 


15.2 






5.7 


8.4 


IRC+ 10216 (S) 


17.8 


10.7 


10.0 


11.6 


5.8 


10.1 


IRC+ 10216(0) 










11.8 


11.8 


IRAS 15082-4808 




8.9 




37.1 


5.0 


10.7 


IRAS15194-5115 


7.2 


11.3 






6.1 


7.7 


lRC+40540 










3.5 


3.5 


Average 


11.2 


11.5 


10.0 


24.4 


6.3 


8.7 



Glassgold (1982) is adopted. The photodissociation radius of a 
parent species (fx, viz., HCN, C2H2, CS) is determined by 



dr 



/x 



•'exp 



and of a daughter (/xd, viz., C2H) by 



^exp ^exp 



dr 



./xd 



(4) 



(5) 



Go is the unshielded photodissociation rate and is the dust 
shielding distance, given by; 



where L is the bolometric luminosity of the star, c the speed 
of light, and {Q) is the flux-averaged momentum transfer effi- 
ciency. The drift velocity also enters in the expression of the 
heating caused by dust-gas collisions. In the self-consistent 
treatment of the energy balance in the circumstellar envelope a 
value of 0.03 was assumed for {Q), whi ch is also retain e d here . 
Photodissocation rates are taken from van Dishoeck ( 1988 ), 
with the assumption that SiS has the same photodissociation 
rate as CS (confirmed to within a factor 2 by the UMIST 
Ratefile^. The calculated shielding distances and rate coeffi- 
cients are shown in Tables ^ and ^, respectively. 



1.4 



3ex Md 
4fldPd 47rvd ' 



(6) 



where Q is the dust absorption efficiency, Md the dust mass- 
loss rate, and Vd the expansion velocity of the dust grains (J ura 
& Morris 1981 ). Since different molecules are generally disso- 
ciated at different wavelengths the adopted value of Q will de- 
pend on the species under study. However, the wavelength de- 
pendence ofgin the region of interest, ~ 1000-3000 A, is weak 
(e.g., Suh 200c ) and a generic value of 1.0 for most species of 
interest here is adopted. However, the shielding distance of CN 
i s taken to be a factor of 1.2 greater than that used for HCN 
( Truong-Bach et al. 1987| ). Introducing the /i-parameter defined 
in Eq. ^ the dust shielding distance may now be expressed as 



hM 

d = 5.27 10^^ cm, 

Vd 



(7) 



where the H2 mass-loss rate M is given in yr ' and Vd in 
km s"'. A gas-to-dust mass ratio of 0.01 was assumed. 

There will g enerally be a drift velocity between the dust 



and the gas (e.g., [Schoier & Olofsson||2001[ ) 



Vd ■ 



^^exp 



i-Vexp<2> 



\ Mc 



(8) 



- iittp: //www. rate99.co.uk 



Table 7. Abundances and sizes of emission regions" 



Molecule IRC+10216 (SEST) IRC+ 10216 (OSO) IRAS 15 194-5 115 IRAS 15082-4808 





n 


r 


e 


/> 




r\ 


r 


e 


/x 


r\ 


r 


e 


/x 


n 


r 


e 


/x 




[cm] 


[cm] 


[X]/[ 


m 


[cm] 


[cm] 


[X]/[H2] 


[cm] 


[cm] 


[X]/[H,] 


[cm] 


[cm] 


[X]/[H2] 


CN(l-O) 


2.6(16) 


7.1 


[16) 


3.4 


(-06) 


2.6(16) 


7.1 


(16) 


2.1 (-06) 


3.5(16) 


9.6 


(16) 


2.1 


'-06) 


2.0(16) 


6.2 


(16) 


3.2 


(-06) 


"CN(l-O) 


2.6(16) 


7.1 


(16) 


7.6 


(-08) 


2.6(16) 


7.1 


(16) 


4.6 (-08) 


3.5(16) 


9.6 


(16) 


3.0 


'-07) 


2.0(16) 


6.2 


(16) 


< 1.3 


(-08) 


CN(2-1) 


2.6(16) 


7.1 


[16) 


9.6 


(-07) 


2.6(16) 


7.1 


(16) 





3.5(16) 


9.6 


(16) 


3.5 


(-07) 


2.0(16) 


6.2 


(16) 


3.8 


(-07) 


CS(2-1) 




4.0 


[16) 


9.9 


(-07) 




4.0 


(16) 


4.6 (-07) 




5.5 


(16) 


2.3 


;-06) 




3.3 


(16) 


2.2 


(-06) 


"CS(2-1) 




4.0 


[16) 


2.2 


(-08) 




4.0 


(16) 







5.5 


(16) 


5.0 


(-07) 




3.3 


(16) 


<3.1 


(-08) 


C3^S(2-1) 




4.0 


[16) 


3.7 


(-08) 




4.0 


(16) 







5.5 


(16) 









3.3 


(16) 







CS(5^) 




4.0 


[16) 


1.2 


(-06) 




4.0 


(16) 







5.5 


(16) 


1.6 


'-06) 




3.3 


(16) 


9.0 


(-07) 


C2S(6,7-5,6) 


1.2(16) 


4.0 


[16) 


<4.6 


(-09) 


1.2(16) 


4.0 


(16) 


<4.5 (-09) 


8.0(15) 


5.5 


(16) 


<2.4 


(-08) 


7.4(15) 


3.3 


(16) 


<3.1 


(-08) 


€28(8,9-7,8) 


1.2(16) 


4.0 


[16) 


<4.1 


(-09) 


1.2(16) 


4.0 


(16) 


<9.5 (-09) 


8.0(15) 


5.5 


(16) 


<3.5 


(-08) 


7.4(15) 


3.3 


(16) 


<4.8 


(-08) 


C3S(15-14) 


1.2(16) 


4.0 


[16) 


< 1.2 


(-08) 


1.2(16) 


4.0 


(16) 


< 1.5 (-08) 


8.0(15) 


5.5 


(16) 


<7.2 


'-08) 


7.4(15) 


3.3 


(16) 


< 1.3 


(-06) 


SiO(2-l) 




2.0 


[16) 


1.3 


(-07) 




2.0 


(16) 


9.6 (-08) 




1.3 


(16) 


1.7 


'-06) 




3.3 


(16) 


7.2 


(-07) 


SiS(5^) 




2.0 


[16) 


1.2 


(-06) 




2.0 


(16) 


8.1 (-07) 




1.3 


(16) 


4.9 


'-06) 




3.3 


(16) 


3.3 


(-06) 


SiS(6-5) 




2.0 


(16) 


9.0 


(-07) 




2.0 


(16) 


9.0 (-07) 




1.3 


(16) 


4.2 


(-06) 




3.3 


(16) 


2.0 


(-06) 


SO(3-2) 


2.6(16) 


7.1 


,16) 


<3.5 


(-09) 


2.6(16) 


7.1 


,16) 





3.5(16) 


9.6 


(16) 







2.0(16) 


6.2 


(16) 







HCN(l-O) 




3.4 


(16) 


1.4 


(-05) 




3.4 


(16) 


1.1 (-05) 




4.6 


(16) 


1.2 


(-05) 




2.7 


(16) 


1.0 


(-05) 


H'3CN(l-0) 




3.4 


(16) 


3.1 


(-07) 




3.4 


(16) 


2.5 (-07) 




4.6 


(16) 


2.0 


(-06) 




2.7 


(16) 


2.9 


(-07) 


HNC(l-O) 


2.4(16) 


8.4 


[16) 


7.2 


(-08) 


2.4(16) 


8.4 


[16) 


3.8 (-08) 


1.6(16) 


5.6 


(16) 


3.3 


(-07) 


1.5(16) 


5.2 


(16) 


1.6 


(-07) 


HN"C(l-0) 


2.4(16) 


8.4 


(16) 


<1.9 


(-09) 


2.4(16) 


8.4 


(16) 





1.6(16) 


5.6 


(16) 


I.l 


(-07) 


1.5(16) 


5.2 


(16) 


<1.2 


(-08) 


SiC2(4,04-3,03) 


2.4(16) 


6.0 


(16) 


1.6 


(-07) 


2.4(16) 


6.0 


(16) 





1.6(16) 


4.0 


(16) 







1.5(16) 


3.7 


(16) 


3.3 


(-07) 


SiC2(4,22-3,21) 


2.4(16) 


6.0 


(16) 


3.2 


(-07) 


2.4(16) 


6.0 


(16) 


— 


1.6(16) 


4.0 


(16) 







1.5(16) 


3.7 


(16) 


<1.8 


(-07) 


SiC2(5,05-4,04) 


2.4(16) 


6.0 


(16) 


1.7 


(-07) 


2.4(16) 


6.0 


(16) 


1.7 (-07) 


1.6(16) 


4.0 


(16) 


1.2 


(-06) 


1.5(16) 


3.7 


(16) 


4.9 


(-07) 


C2H(l-0) 


2.3(16) 


5.6 


(16) 


2.8 


(-06) 


2.3(16) 


5.6 


(16) 


2.4 (-06) 


1.6(16) 


4.8 


(16) 


1.5 


(-05) 


1.7(16) 


4.8 


(16) 


8.3 


(-06) 


C3H(9/2-7/2) 


2.4(16) 


8.4 


(16) 


5.5 


(-08) 


2.4(16) 


8.4 


(16) 


5.9 (-08) 


1.6(16) 


5.6 


(16) 


9.2 


'-08) 


1.5(16) 


5.2 


(16) 


< 1.7 


(-08) 


C3N(9-8) 


2.4(16) 


8.4 


(16) 


3.0 


(-07) 


2.4(16) 


8.4 


(16) 





1.6(16) 


5.6 


(16) 


<4.9 


(-08) 


1.5(16) 


5.2 


(16) 


<1.0 


(-07) 


C3N(11-10) 


2.4(16) 


8.4 


(16) 


5.9 


(-07) 


2.4(16) 


8.4 


(16) 


7.6 (-07) 


1.6(16) 


5.6 


(16) 


9.2 


(-07) 


1.5(16) 


5.2 


(16) 


9.5 


(-07) 


C4H(9-8) 


2.4(16) 


8.4 


(16) 


3.7 


(-06) 


2.4(16) 


8.4 


(16) 





1.6(16) 


5.6 


(16) 







1.5(16) 


5.2 


(16) 


<6.6 


(-07) 


C4H(10-9) 


2.4(16) 


8.4 


(16) 


2.7 


(-06) 


2.4(16) 


8.4 


(16) 


3.2 (-06) 


1.6(16) 


5.6 


(16) 


2.8 


(-05) 


1.5(16) 


5.2 


(16) 


6.9 


(-06) 


C3H2(2,12-1,01) 


2.4(16) 


8.4 


(16) 


3.2 


(-08) 


2.4(16) 


8.4 


(16) 


— 


1.6(16) 


5.6 


(16) 


6.0 


'-07) 


1.5(16) 


5.2 


(16) 


<9.3 


(-08) 


HC3N(10-9) 


1.2(16) 


6.0 


(16) 


1.1 


(-06) 


1.2(16) 


6.0 


(16) 


4.4 (-07) 


8.0(15) 


4.0 


(16) 


1.9 


'-06) 


7.4(15) 


3.7 


(16) 


2.3 


(-06) 


HCC"CN(10-9) 


1.2(16) 


6.0 


(16) 


2.4 


(-08) 


1.2(16) 


6.0 


(16) 


< 2.8 (-09) 


8.0(15) 


4.0 


(16) 


3.9 


'-07) 


7.4(15) 


3.7 


(16) 


<2.7 


(-08) 


HC"CCN(10-9) 


1.2(16) 


6.0 


(16) 


2.4 


(-08) 


1.2(16) 


6.0 


(16) 


< 2.8 (-09) 


8.0(15) 


4.0 


(16) 


3.9 


(-07) 


7.4(15) 


3.7 


(16) 


<2.7 


(-08) 


HC3N(11-10) 


1.2(16) 


6.0 


(16) 


7.5 


(-07) 


1.2(16) 


6.0 


(16) 





8.0(15) 


4.0 


(16) 






7.4(15) 


3.7 


(16) 






HC3N( 12-11) 


1.2(16) 


6.0 


(16) 


1.9 


(-06) 


1.2(16) 


6.0 


(16) 


— 


8.0(15) 


4.0 


(16) 


2.2 


'-06) 


7.4(15) 


3.7 


(16) 


2.1 


(-06) 


HCC"CN(12-11) 


1.2(16) 


6.0 


(16) 


4.1 


(-08) 


1.2(16) 


6.0 


(16) 


1.5 (-08) 


8.0(15) 


4.0 


(16) 


2.8 


(-07) 


7.4(15) 


3.7 


(16) 


<2.0 


(-08) 


HC"CCN(12-11) 


1.2(16) 


6.0 


(16) 


4.1 


(-08) 


1.2(16) 


6.0 


(16) 


1.5 (-08) 


8.0(15) 


4.0 


(16) 


2.8 


(-07) 


7.4(15) 


3.7 


(16) 


<2.0 


(-08) 


CH3CN(6(1)-5(1)) 


1.2(16) 


6.0 


(16) 


1.2 


(-08) 


1.2(16) 


6.0 


(16) 


< 2.0 (-09) 


8.0(15) 


4.0 


(16) 


<2.5 


'-08) 


7.4(15) 


3.7 


(16) 


<4.3 


(-08) 


CH3CN( 12(0)- 11(0)) 


1.2(16) 


6.0 


,16) 


< 1.2 


(-07) 


1.2(16) 


6.0 


,16) 




8.0(15) 


4.0 


(16) 


<2.4 


(-06) 


7.4(15) 


3.7 


(16) 


<1.5 


(-06) 


HC5N(32-31) 


1.2(16) 


6.0 


(16) 


3.9 


(-06) 


1.2(16) 


6.0 


(16) 




8.0(15) 


4.0 


(16) 


<4.5 


(-06) 


7.4(15) 


3.7 


(16) 


<3.8 


(-06) 


HC5N(34-33) 


1.2(16) 


6.0 


(16) 


4.2 


(-06) 


1.2(16) 


6.0 


(16) 


< 7.6 (-07) 


8.0(15) 


4.0 


(16) 






7.4(15) 


3.7 


(16) 


<7.4 


(-06) 


HC5N(35-34) 


1.2(16) 


6.0 


[16) 


1.3 


(-05) 


1.2(16) 


6.0 


(16) 




8.0(15) 


4.0 


(16) 






7.4(15) 


3.7 


(16) 


<1.1 


(-05) 



(a) In this table x (y) represents xxW. 

(b) Bold face indicates an abundance calculated by scaling the abundance of the '^C isotope of the same species by the modelled '^CO/'^CO-ratio (see Sect. 5.5). 



Table 7. (cont.) Abundances and sizes of emission regions' 



Molecule IRAS07454-7112 CIT6 AFGL3068 IRC+40540 





n 


r 


e 


fx 


n 


re 


fx 


n 


r 


e 


fx 


n 


r 


e 


fx 




cm 


cm 


[X]/[ 


H2] 


cm 


cm 


[X]/[H2] 


cm 


cm 


[X]/[H,] 


cm 


cm 


[X]/[H2] 


CN(l-O) 


1.3(16) 


4.1 


(16) 


5.3 


[-06) 


1.6(16) 


5.1(16) 


2.0 (-05) 


3.2(16) 


8.3 


(16) 


4.6 (-07) 


3.2(16) 


8.3 


(16) 


1.0 (-06) 


"CN(l-O) 


1.3(16) 


4.1 


(16) 


4.8 


(-07) 


1.6(16) 


5.1(16) 


5.8 (-07) 


3.2(16) 


8.3 


(16) 


< 2.6 (-08) 


3.2(16) 


8.3 


(16) 


< 9.9 (-09) 


CN(2-1) 


1.3(16) 


4.1 


(16) 


1.6 


[-06) 


1.6(16) 


5.1(16) 





3.2(16) 


8.3 


(16) 





3.2(16) 


8.3 


(16) 





CS(2-1) 




2.2 


(16) 


1.6 


[-06) 




2.7(16) 


2.5 (-06) 




4.8 


(16) 


3.7 (-07) 




4.9 


(16) 


5.4 (-07) 


"CS(2-1) 




2.2 


(16) 


<2.9 


[-08) 




2.7(16) 







4.8 


(16) 







4.9 


(16) 





C^^S (2-1) 




2.2 


(16) 









2.7(16) 







4.8 


(16) 







4.9 


(16) 





CS(5^) 




2.2 


(16) 


1.6 


[-06) 




2.7(16) 


— 




4.8 


(16) 







4.9 


(16) 





€28(6,7-5,6) 


5.6(15) 


2.2 


(16) 


<9.1 


[-08) 


4.3(15) 


2.7(16) 


< 1.0 (-07) 


1.6(16) 


4.8 


(16) 


< 8.1 (-08) 


1.6(16) 


4.9 


(16) 


< 3.0 (-08) 


€28(8,9-7,8) 


5.6(15) 


2.2 


(16) 


<1.1 


[-07) 


4.3(15) 


2.7(16) 


< 2.0 (-07) 


1.6(16) 


4.8 


(16) 


< 8.0 (-08) 


1.6(16) 


4.9 


(16) 


< 5.6 (-08) 


€38(15-14) 


5.6(15) 


2.2 


(16) 


<3.6 


[-07) 


4.3(15) 


2.7(16) 


< 2.5 (-07) 


1.6(16) 


4.8 


(16) 


< 3.5 (-07) 


1.6(16) 


4.9 


(16) 


< 8.9 (-08) 


SiO(2-l) 




9.4 


(15) 


4.4 


[-07) 




7.2(15) 


1.0 (-06) 




2.6 


(16) 


<4.7 (-08) 




2.6 


(16) 


5.6 (-08) 


818(5^) 




9.4 


(15) 


4.8 


[-06) 




7.2(15) 


1.0 (-06) 




2.6 


(16) 


3.3 (-07) 




2.6 


(16) 


1.4 (-06) 


818(6-5) 




9.4 


(15) 


3.4 


[-06) 




7.2(15) 


5.2 (-06) 




2.6 


(16) 


1.0 (-06) 




2.6 


(16) 


4.9 (-07) 


HCN(l-O) 




1.8 


(16) 


7.8 


(-06) 




2.2(16) 


1.1 (-05) 




4.2 


(16) 


6.3 (-06) 




4.2 


(16) 


6.5 (-06) 


H"CN(l-0) 




1.8 


(16) 


5.6 


(-07) 




2.2(16) 


3.0 (-07) 




4.2 


(16) 


2.1 (-07) 




4.2 


(16) 


1.3 (-07) 


HNC(l-O) 


1.1(16) 


4.0 


(16) 


1.0 


(-07) 


8.6(15) 


3.0(16) 


2.3 (-07) 


3.1(16) 


1.1 


(17) 


3.0 (-08) 


3.1(16) 


1.1 


(17) 


2.2 (-08) 


HN"C(l-0) 


1.1(16) 


4.0 


(16) 


<1.8 


(-08) 


8.6(15) 


3.0(16) 


< 2.1 (-08) 


3.1(16) 


1.1 


(17) 


< 1.4 (-08) 


3.1(16) 


1.1 


(17) 


< 4.4 (-09) 


8iC2(4,04-3,03) 


1.1(16) 


2.8 


(16) 


2.3 


(-07) 


8.6(15) 


2.2(16) 





3.1(16) 


7.8 


(16) 





3.1(16) 


7.8 


(16) 





8iC2(4,22-3,21) 


1.1(16) 


2.8 


(16) 


<2.7 


(-07) 


8.6(15) 


2.2(16) 





3.1(16) 


7.8 


(16) 





3.1(16) 


7.8 


(16) 





8iC2(5,05^,04) 


1.1(16) 


2.8 


(16) 


<3.2 


(-07) 


8.6(15) 


2.2(16) 


3.1 (-06) 


3.1(16) 


7.8 


(16) 


< 7.5 (-08) 


3.1(16) 


7.8 


(16) 


< 5.3 (-08) 


C2H(l-0) 


1.1(16) 


3.2 


(16) 


<4.3 


(-07) 


1.4(16) 


4.1(16) 


6.9 (-06) 


2.8(16) 


6.7 


(16) 


5.7 (-06) 


2.7(16) 


6.7 


(16) 


< 1.0 (-07) 


C3H(9/2-7/2) 


1.1(16) 


4.0 


(16) 


<3.1 


(-08) 


8.6(15) 


3.0(16) 


< 7.4 (-08) 


3.1(16) 


1.1 


(17) 


< 1.4 (-08) 


3.1(16) 


1.1 


(17) 


< 1.2 (-08) 


C3N(9-8) 


1.1(16) 


4.0 


(16) 


<1.6 


(-07) 


8.6(15) 


3.0(16) 





3.1(16) 


1.1 


(17) 





3.1(16) 


1.1 


(17) 





C3N(11-10) 


1.1(16) 


4.0 


(16) 


7.3 


(-07) 


8.6(15) 


3.0(16) 


2.6 (-06) 


3.1(16) 


1.1 


(17) 


5.5 (-07) 


3.1(16) 


1.1 


(17) 


1.1 (-07) 


C4H(9-8) 


1.1(16) 


4.0 


(16) 







8.6(15) 


3.0(16) 





3.1(16) 


1.1 


(17) 





3.1(16) 


1.1 


(17) 





C4H(10-9) 


1.1(16) 


4.0 


(16) 


<8.8 


(-07) 


8.6(15) 


3.0(16) 


< 1.7 (-06) 


3.1(16) 


1.1 


(17) 





3.1(16) 


1.1 


(17) 





C3H2(2,12-1,01) 


1.1(16) 


4.0 


(16) 







8.6(15) 


3.0(16) 





3.1(16) 


1.1 


(17) 





3.1(16) 


1.1 


(17) 





HC3N(10-9) 


5.6(15) 


2.8 


(16) 


1.7 


(-06) 


4.3(15) 


2.2(16) 


2.4 (-06) 


1.6(16) 


7.8 


(16) 


5.0 (-07) 


1.6(16) 


7.8 


(16) 


1.5 (-06) 


HCC"CN(10-9) 


5.6(15) 


2.8 


(16) 


1.5 


(-07) 


4.3(15) 


2.2(16) 


< 4.9 (-08) 


1.6(16) 


7.8 


(16) 


< 1.5 (-08) 


1.6(16) 


7.8 


(16) 


3.1 (-08) 


HC"CCN(10-9) 


5.6(15) 


2.8 


(16) 


1.5 


[-07) 


4.3(15) 


2.2(16) 


< 4.9 (-08) 


1.6(16) 


7.8 


(16) 


< 1.5 (-08) 


1.6(16) 


7.8 


(16) 


3.1 (-08) 


HC3N(11-10) 


5.6(15) 


2.8 


(16) 




— 


4.3(15) 


2.2(16) 





1.6(16) 


7.8 


(16) 





1.6(16) 


7.8 


(16) 





HC3N(12-11) 


5.6(15) 


2.8 


(16) 


2.4 


(-06) 


4.3(15) 


2.2(16) 


— 


1.6(16) 


7.8 


(16) 


— 


1.6(16) 


7.8 


(16) 


— 


HCC"CN(12-11) 


5.6(15) 


2.8 


(16) 


1.3 


(-07) 


4.3(15) 


2.2(16) 


2.7 (-07) 


1.6(16) 


7.8 


(16) 


< 2.9 (-08) 


1.6(16) 


7.8 


(16) 


< 1.1 (-08) 


HC"CCN(12-11) 


5.6(15) 


2.8 


(16) 


1.3 


(-07) 


4.3(15) 


2.2(16) 


2.7 (-07) 


1.6(16) 


7.8 


(16) 


< 2.9 (-08) 


1.6(16) 


7.8 


(16) 


< 1.1 (-08) 


CH3CN(6(1)-5(1)) 


5.6(15) 


2.8 


(16) 


<2.4 


(-07) 


4.3(15) 


2.2(16) 


< 1.3 (-07) 


1.6(16) 


7.8 


(16) 


< 1.2 (-07) 


1.6(16) 


7.8 


(16) 


< 6.4 (-08) 


CH3CN(12(0)-11(0)) 


5.6(15) 


2.8 


(16) 


<4.6 


(-06) 


4.3(15) 


2.2(16) 




1.6(16) 


7.8 


(16) 




1.6(16) 


7.8 


(16) 




HC5N(32-31) 


5.6(15) 


2.8 


(16) 






4.3(15) 


2.2(16) 




1.6(16) 


7.8 


(16) 




1.6(16) 


7.8 


(16) 




HC5N(34-33) 


5.6(15) 


2.8 


(16) 


9.2 


(-06) 


4.3(15) 


2.2(16) 


< 1.3 (-05) 


1.6(16) 


7.8 


(16) 


< 4.1 (-06) 


1.6(16) 


7.8 


(16) 


< 4.5 (-06) 


HC5N(35-34) 


5.6(15) 


2.8 


(16) 


<8.8 


(-06) 


4.3(15) 


2.2(16) 




1.6(16) 


7.8 


(16) 




1.6(16) 


7.8 


(16) 





(a) In this table x (y) represents xxW. 

(b) Bold face indicates an abundance calculated by scaling the abundance of the '^C isotope of the same species by the modelled '^CO/'^CO-ratio (see Sect. 5.5). 
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The photodissociation radii of species formed in the pho- 
tosphere were chosen to be the e-folding radii of the initial 
abundances, i.e. /x(?"e) = /x(^*)/e, where R^, is the stellar ra- 
dius. For the photodissociation products inner and outer radii 
are chosen to be where the abundance has dropped to 1/e of its 
peak value. The envelope sizes calculated by the simple pho- 
todissociation model agree relatively well with observed val- 
ues in IRC -1-10216 and other objects, except in the case of SiS 
(see next paragraph). Lindqvist et al. (2000) used a detailed ra- 
diative transfer method and observed brightness distributions 
to derive envelope sizes of HCN and CN. For IRC -1-102 16 and 
IRC+40540 they found envelope sizes of -4x10"" cm for HCN 
and ~5x 10'^ cm for CN, in excellent agreement with the values 
derived here from the photodissociation model. In the case of 
CIT 6 the observed HCN envelope size is larger than that cal- 
culated by a factor of 2. Envelope sizes and derived abundances 
are given in Table ^ 

SiO and SiS are parent species with a radial extent depend- 
ing on photodissociation. The calculation of the SiS photodis- 
sociation radius in IRC -1-10216, however, is not consistent with 



the radius observed by Bieging & Tafalh (1993). Therefore the 
observed radius for SiS in lRC-i-10216 was used in the abun- 
dance calculations, and its radius was scaled for the other ob- 
jects in the sample with the factor (M/vex£)*/(M/vexp)io2i6 in 



the same way as described in section p.4.2| . The outer radius 
calculated via the photodissociation model is actually slightly 
larger than that observed. This gives credence to the idea that 
SiS freezes out onto dust grains. SiO is treated in the same way 
as SiS since it too is likely to freeze out onto grains. Hence the 
same envelope size as SiS was adopted for this species. 

5.4.2. Inner and outer radii for species of circumstellar 
chemistry origin 

For those molecules that have their origin in a circumstellar 
chemistry a different approach is used. The sizes of the emis- 
sion regions of HNC and HC3N have been determined interfer- 
ometr ically for IRC-H10216 using observation s by Bieging & 
Rieu ( |1988| ), and of SiCi by j^^ucasetalj ( |1995| ). Inner and outer 
radii for HNC, HC3N and SiCa are taken from observations 
of IRCh-10216, and scaled by a factor (M/vexp)*/(M/vexp)i{)2i6 
for the remaining six sample carbon stars (since molecular dis- 
tributions have not been mapped in all the sample stars; see 
Table ^ for values of this ratio). This factor is proportional to 
the density ratio at a given radius, and also (to the first order) 
to the ratio of the envelope shielding distances. Chemical mod- 
elUng shows that the species C3H, C4H, C3H2, C3N, CH3CN, 
HC5N, SO, C2S and C3S are formed via chemical reactions in 



the outer envelope (e.g., Millar & Herbst 1994). Furthermore, 
it shows that the species C3H, C4H, C3H2 and C3N have a simi- 
lar radial distribution to HNC, and that the species CH3CN and 
HC5N are similarly distributed to HC3N. Hence corresponding 
radii for these species are assumed in the calculations. SO is as- 



sumed to have a similar distribution to CN (c.f., Nejad & Millar 



equal to that of HC3N). All isotopes (i.e. species involving '^C 
and •''*S) are further assumed to have the same distributions as 
the main isotope. 

5.5. Molecular abundances 

5.5.1 . Uncertainties in tine abundance estimates 

The assumption of optically thin emission has a systematic 
effect on the derived abundances, in that the true abundance 
will be higher if there are opacity effects. As can be seen from 
Table there can be a factor of 2-3 in error from this assump- 
tion. 

The accuracy of the calculated abundances is dependent on 
various assumptions. Due to the intrinsic difficulties in estimat- 
ing distances to the objects included in the sample, the typi- 
cal uncertainty in the adopted distance is a factor of ~2. This 
influences the mass-loss rates derived in the radiative transfer 
modelling such that the adopted mass-loss rate will scale as 
D'"^. The calculation of the limiting radii of a certain molec- 
ular distribution in the envelope is also affected by inaccura- 
cies in the distance estimate. Where radii have been calculated 
by scaling observed radii in lRC-t-10216, there comes an er- 
ror which scales as D. The results from the photodissociation 
model will have a lesser dependence, with D coming into the 
expression for the shielding distance, via M. Hence the abun- 
dance, which is trivially derived from Eq. ^, will vary as ap- 
proximately D"'"", giving a factor of 2, possibly, in error 

In the present analysis an excitation temperature of lOK 
was assumed for all transitions. The error in abundance esti- 
mate due to this assumption will depend on the excitation tem- 
perature of a particular transition and the assumption of LTE, 
and is estimated to be not more than a factor of ~2. 

Overall, it seems like an error of a factor of 5 is reasonable 
to expect in the abundances presented here, although it should 
be borne in mind that this could possibly rise to an order of 
magnitude. In the comparison of abundance ratios in Table ^ 
any difference less than a factor of 5 is treated as insignificant 
with respect to error margins. 

5.5.2. Comparison witli previously published 
abundances 



In comparison with jVyman et al. ( 1993 ), derived abundance 
ratios (Table H) have increased in favour of IRAS 15194- 
5115 by up to a factor of approximately 4. Individual abun- 
dances generally show a factor ~2 increase for those in IRAS 
15914-5115, and a factor ~2 decrease for IRC-fl0216, over 



1988), and C2S and C3S are assumed to follow the CS distribu- 



tion (c.f., Millar et al. 2001), although these latter two species 



are not parent molecules (and hence are given an inner radius 



Nyman et alj (|1993|). In |Nyman et al.| ( |1993[ ), the distances 
adopted for IRC -Hi 02 16 and IRAS 15194-5115 were approx- 
imately twice as large as those here; however this has little ef- 
fect on calculated abundances since the distance ratio is more 
or less unchanged. However, the calculated mass-loss rate of 
IRAS 15 194-5 1 15 in the previous paper was larger than that of 
IRC -1-10216 by a factor of 2.5, and here the recalculated mass- 
loss rates are the same for these two objects. This would tend 
to increase the abundance ratios quoted by a similar factor The 
photodissociation radius depends on the mass-loss rate and the 
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Table 8. Molecular isotope abundance ratios 





1RC+10216(S) 


lRC+10216 (0) 


IRAS 15 194 


1RAS07454 


C1T6 


lRC+40540 


HC3N/HC'''CCN(10-9) 


32.6 




4.9 


11.1 




17.9 


HC3N/HC'3CCN(12-11) 


22.1 




7.9 


17.8 






CS/'3CS(2-1) 


22.7 




4.6 








CN/"CN(l-0) 


25.1 


46.6 


6.7 


11.0 


12.4 




Average 


25.6 


46.6 


6.0 


13.3 


12.4 


17.9 


Modelling of CO (Table 


45.0 


45.0 


6.0 


14.0 


35.0 


50.0 



dust parameters through the dust shielding distance. Compared 
to Nymanetal. (1993) this paper uses different mass-loss rates, 
a different gas-to-dust ratio (0.01 compared to 0.005), and the 
dust parameters are also slightly different due to the determina- 
tion of the /i-factor in the radiative transfer analysis. The scale 
factor that determines the inner and outer radii of the species 
with a origin in circumstellar chemistry depends on the ratio 
of the mass-loss rates. Thus the difference in relative mass loss 
rates explains the factor of 4 difference in relative abundances 
between IRAS 15914-51 15 and I RC+ 102 16 derived in this pa- 
per compared to those derived in Nyman et al. (1993). 



Table 10. Fractional abundances for IRC -1-102 16 for species of 
photospheric origin. 



As discussed earlier in this paper the outer radii of SiS 
and SiO are not determined through the photodissociation ra- 
dius but scaled from their observed radius in IRC -(-10216. In 
this way the calculated abundances of SiO and SiS have in- 
creased in by a factor of 4-5 for IRAS 15194-51 15 compared 
to Nyman et al. (1993 ). For IRC-fl0216 these species have the 
same abundance as calculated previously, and they are in rea- 
sonable agreement with those reported elsewhere in the litera- 

: 5.6 10-^ and /(SiS) 



ture. 



. Bujarrabalet al.| ( |1994| ) give /(SiO) ■■ 
- 3.9 10"*', which are ~4 times greater, using an outer radius 
half that quoted in this paper and a distance of 200 pc. Better 



agreement is seen for the other species which Bujarrabal et al 



(1994) detect, with the exception of HNC, which is a factor 
~6 in disagreement. [Cernicharo, Guelin, & Kahane ( 200C ) also 
calculate abundances in IRC+10216, and are within a factor 5 
of those here. 



Generally, it seems that IRC+10216 has lower abundances 
than IRAS 15194-5115, IRAS 15082-4808, IRAS 07454- 
7112 and CIT 6. However, it is very similar, physically and 
chemically, to AFGL 3068 and IRC-H40540. 

The three northern sources observed at OSO have also 
been studied in the literature, and abundances derived. The 



Bujarrabal et al, (1994) paper includes data relating to CIT 
6, AFGL 3068 and IRCh-40540. All abundances calculated by 
Bujarrabal et al. in these three sources are greater than those de- 
rived here. Distances and mass-loss rates are reasonably com- 
parable between this paper and that. Generally, this means that 
the calculated abundances in this paper are lower by a factor of 
~2 in CIT 6, a factor of ~5 in AFGL 3068, and a factor of ~3 
in IRC -1-40540 compared to those derived in Bujarrabal et al. 
(1994). 



Species 


Observed" 


Model'' 






TE Shock 


CO 


1.0(-3)'^ 


9.8(-4) 1.0(-3) 


HCN 


1.3(-5) 


5.1(-5) 3.1(-6) 


CS 


9.9(-7) 


1.3(-5) 7.1(-7) 


SiS 


9.5(-7) 


1.0(-5) 2.5(-5) 


SiO 


1.1 (-7) 


1.9(-8) 2.7(-7) 



(a) From Table ^ 



(b) From chemical modelling by Willacy & Cherchnefi ( 1998) 



(c) Assumed value that determines molecular hydrogen densities. 

6. Discussion 

6.1. Chemistry 

6.1 .1 . Presence of shocks in the inner wind? 

When thermal equilibrium (TE) prevails the molecular content 
of a gas can be readily calculated from its elemental chemi- 
cal composition. This is the case in the stellar photosphere and 
near the inner boundary of the envelope of an AGB-star, where 
the gas density and temperature are high. The variable nature 
of AGB-stars induces pulsation-driven shocks that propagate 
outwards and suppress TE from the point of shock formation. 
Non-equilibrium chemical m odelling has been performed by 
Willacy & Cherchnefi} ( [l99l for IRC-H10216), and shows that 
shocks can strongly alter the chemical abundances in the inner 
regions of the CSE from their TE values. In particular, SiO is 
strongly enhanced whereas HCN and CS are destroyed. Other 
species, e.g., CO and SiS are relatively unaffected. 

In Table the abundances obtained by Willacy & 



Cherchneff (1998) for TE and non-equilibrium chemical mod- 
elling (shock strength of 1 1.7 kms"') are compared to the val- 
ues obtained in the analysis. Given the uncertanties, about a 
factor of five, the abundances obtained from the observations 
clearly favour a scenario in which a shock has passed through 
the inner (< 5 R*) parts of the wind. The SiO and SiS abun- 
dances which are derived are lower limits since they could be 
significantly depleted in the outer envelope due to freeze-out 
onto dust grains. The average abundances derived for HCN, 
CS, and SiO in the sample are 1 x 10"^ 1.5 x 10"^ 7 x 10"^ re- 
spectively. The SiO abundance shows the largest spread among 
the sources reflecting its sensitivity to the shock strength and 
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Molecule 


IRC+10216 
SEST OSO 


IRAS 15 194 


IRAS 15082 


IRAS07454 


CIT6 


AFGL3068 


IRC+40540 


CN(l-O) 


1 n 


i. i 


1.1 


1.7 


2.8 


3.8 


0.2 


0.5 


"CN(l-O) 


i.U 


U.O 


4.0 




6.3 


7.6 






CN(2-1) 


I.U 




0.4 


0.4 


1.7 








CS(2-1) 


I.U 


u.y 


4.6 


4.3 


3.2 


5.1 


0.7 


1.1 


'''CS(2-1) 


1 n 
i .U 




22.7 












CS(5-4) 


1 n 
I.U 




1.4 


0.8 


1.4 








SiO(2-l) 


1 n 
i.U 


U. / 


12.8 


5.4 


3.3 


7.8 




0.4 


SiS(5-4) 


1 n 

i.U 


n 7 
u. / 


4.0 


2.7 


3.9 


0.8 


0.3 


1.1 


SiS(6-5) 


i.U 


1 n 
I.U 


4.7 


2.2 


3.8 


5.7 


1.1 


0.5 


HCN 


i.U 


n Q 
U.5 


0.9 


0.7 


0.6 


0.8 


0.5 


0.5 


H"CN(l-0) 


1 n 
i.U 


n c 
U.5 


6.4 


0.9 


1.8 


1.0 


0.7 


0.4 


HNC(l-O) 


1 n 
i.U 


U.3 


4.6 


2.2 


1.4 


3.2 


0.4 


0.3 


SiC2(4,04-3,03) 


1 n 
i.U 






2.0 


1.5 








SiC2(5,05-4,04) 


1 n 
i.U 


I.U 


6.9 


2.8 




17.7 






C2H(l-0) 


1 n 
i.U 


U.V 


5.6 


3.0 




2.5 


2.0 




CiH(9/2-7/2) 


1.0 


1.1 


1.7 












C3N(11-10) 


1.0 


1.3 


1.5 


1.6 


1.2 


4.4 


0.9 


0.2 


C4H(10-9) 


1.0 


1.2 


10.4 


2.6 










C3H2(2,12-1,01) 


1.0 




18.6 












HC3N(10-9) 


1.0 


0.6 


2.4 


3.0 


2.2 


3.1 


0.6 


0.7 


H'"'C3N(10-9)'' 


1.0 




16.3 




6.4 






1.3 


HC3N(12-11) 


1.0 




2.5 


2.3 


2.6 








H<"'C3N(12-11)'' 


1.0 


0.4 


6.9 




3.3 


6.6 






HC5N(34-33) 


1.0 








2.2 









Bold face signifies a factor of more than 5. 

(a) signifies blend of HCC'^CN and HC^CCN. 



possible variation in the C/O-ratio among the sample sources. 
CS might, however, not be particularly well suited as a probe of 
shocked chemistry. Olofsson et al, (1993c) found, when mod- 
elling a large sample of optically bright carbon stars, that in 
their photospheric LTE models the CS abundance varied con- 
siderably with adopted stellar temperature. 

6.1 .2. Photochemistry in the outer envelope 

Many photochemical models have been developed for the outer 
envelope of IRC+10216 (e.g. [Millar & Herbst||l994t Doty & 



are similar When comparing the present sample of carbon stars 
to the sample of Olofsson et al. (1993a), which on the whole 
tend to have a low mass-loss rate (~10"^Moyr"'), the agree- 
ment in derived abundances for the star in common, CIT 6, is 
very good. There is less than a factor 3 difference in the calcu- 
lated abundances of HCN, CN and CS. However, in general, 
the sample of low mass-loss rate stars has calculated abun- 
dances which are systematically an order of magnitude greater 
than those calculated here. It must be noted that plofsson et al 



(1993a) use an excitation temperature of 20 K, twice that used 



in this analysis. Moreover, Schoier &l Olofsson (2001) show 



Leung 1998; Millar et al. 2000). The physical conditions in the 
outer parts of the wind (> lOOR*) allow for the penetration 
of ambient ultraviolet radiation that induces a photochemistry. 
In Table [IT] the derived column densities for IRC+10216 for a 
number of species produced in the envelope are compared to 



that the mass-l o ss rate s in the low mass-loss rate objects in 
[olofsson et al. (1993a) are underestimated by about a factor 
of 5 on average. This would explain the apparent discrepancy 
in calculated abundances between the high mass-loss rate ob- 



jects here and the lower mass-loss rate objects in Olofsson et al 



those from the photochemical model of [Millar et al.[ (pOOOj). In (li993a|). Hence there seems not to be a marked difference in the 



general the abundances agree well, given the uncertainties. 

It is remarkable that the abundances generally show rela- 
tively little variation within the sample. Most apparent is the 
over-abundance of Si-bearing molecules in CIT 6 (Sect. [6.3[ ). 
Some of the abundances of IRAS 15194-5115 also stand out 
and will be separately discussed in Sect. [6.2[ There are no ap- 
parent trends with the stellar or circumstellar parameters. This 
would suggest that the physical structure in these sources in- 
deed is much the same and that the initial atomic abundances 



molecular composition of high and low mass-loss carbon stars. 

The large spread in the abundance of isotopomers contain- 
ing '■'C follows from the varying '^C/'-'C-ratio among the sam- 
ple sources. The ^^C/'-'C-ratio is related to the nucleosynthesis 
rather than the chemistry and reflects the evolutionary status 
of these stars. However, chemical fractionation may affect this 
ratio in certain molecules, in particular CO. 

CN/HCN and HNC/HCN ratios can also be used to esti- 
mate the evolutionary status of carbon stars. CN is produced 
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via the photodissociation of HCN by ultraviolet radiation, and 
as stellar radiation increases with evolution from AGB star to 
PPN to PN, so the CN/HCN ratio will increase (e.g., Bachiller 
et al. |1997a| ; |Cox et al.||1992D . HNC, formed from the dissocia- 
tive recombination of HCNH+, behaves in a similar way. In this 
sample there is a rather large spread in the CN/HCN ratio, from 
0.07-0.68, and a value of 1.82 for CIT 6 (see Sect. O). This 



large spread is in contrast to Olofsson et al. (1993a), for ex- 
ample, who found a very narrow range (CN/HCN ~0. 65-0. 70) 
in low mass-loss stars. However, there is excellent agreement 
with the results of Lindqvist et al. (2000). For the carbon stars 



IRC+10216, CIT 6 and IRC-t-40540 they derive ratios of 0.16, 
1.5 and 0.17, respectively, in comparison with the 0.22, 1.82 
and 0.15 derived h ere. A CN/HCN ratio of ~0.5 is typical in 
C-rich AGB stars ( Bachiller et al. 1997b ), increasing to ~5 in 
PPNe. In fact, a ratio of 0.6-0.7 is predicted by the photodis- 
sociation model, with only a weak depe ndence on mass-loss 
rate (see Fig. 8 of [Lindqvist et al.||200C| ). The HNC/HCN ra- 
tio seems to be split into two ranges in the present sample of 
carbon stars. IRC-H10216, AFGL 3068 and IRC-H40540 have 
HNC/HCN ratios of <0.005, whilst the remaining stars have 
ratios of 0.01-0.03. The value derived for IRC +10216 is in 
agreement with that quoted in Cox et al. (1992). This seems 



to indicate that the sample stars are not well evolved, since a 
HNC/HCN ratio of ~1 is expected in PPNe (e.g., |Cox et al 



1992D . Having said this, the HNC/HCN ratio does rapidly be- 
come of the order 1 , as can be seen in models of PPN chemistry 
((Woods et al.||2002t). 



Generally, it seems that to use the term "carbon chemistry" 
to refer to a paradigm of chemistry in C-rich evolved stars 
is reasonable. Of the sample stars here, given the variety of 
molecular species, there is very little difference in molecular 
abundances, save for two slightly curious sources, as detailed 
in the following subsections. 



Table 11. Radial column densities (cm ^) for species of cir- 
cumstellar origin towards IRCh-10216. 



Species 


Observed" 


MHB" 


Obs./MHB 


CN 


8.3(14) 


1.0(15) 


0.8 


HNC 


2.0(13) 


8.4(13) 


0.2 


C2H 


8.9(14) 


5.7(15) 


0.2 


C3H 


2.1(13) 


1.4(14) 


0.2 


C3N 


2.0(14) 


3.2(14) 


0.6 


C4H 


9.3(14) 


1.0(15) 


0.9 


C3H2 


1.2(13) 


2.1(13) 


0.6 


HC3N 


9.1(14) 


1.8(15) 


0.5 


CH3CN 


9.9(12) 


3.4(12) 


2.9 


HC5N 


5.8(15) 


7.1(14) 


8.2 



(a) Calculated from Tables | & ^ 

(b) From chemical modelling by Millar et al. ( 200(]| ). 



6.3. CIT 6 

CIT 6 is another object outstanding in the sample. It has a 
CN/HCN ratio of ~1.8, which suggests an advanced evolution- 
ary status, but a low HNC/HCN ratio (~0.02), which suggests 
the contrary. Certainly the idea that CIT 6 is well on its way to 
becoming a PPN has been put forward before (e.g., Trammell 



et al. |1994| ; [Monnier et al. |2000| ; [Zacs et alj|2001[ ). 

The modelled '^CO/'^CO ratio in this source agrees well 



with that carried out previously (Groenewegen et al. 1996 



[Schoier & Olofsson 2000). This ratio, however, does not 



agree with C/ C ratios derived from observations of other 
molecules and their '^C isotopes, both in this paper (Table 



and elsewhere (Kahane et al. 1992; Groenewegen et al. 1996| 

A further point worth note is the comparative over- 
abundance of Si-bearing species which possibly indicates a less 
efficient freeze-out onto dust grains in this particular source. 



6.2. IRAS 15194-5115 



7. Conclusions 



Of the derived abundances those obtained for IRAS 15194- 
5115 stand out the most. In particular, the SiO and C4H abun- 
dances, in addition to the isotopomers containing '^C, ap- 
pear significantly enhanced towards this source. C3H2 also ap- 
pears to be greatly enhanced in this source, but, however, this 
molecule is only observed in one other star, IRC+10216. The 
i^CO/'^CO-ratio of 6 derived for IRAS 15194-5115 is signif- 
icantly lower than that of the others and that which is com- 
monly derived for carbon stars. This value is certain, with the 
modelling of the CO emission being supported by intensity ra- 
tios for another four species, which agree to ±30 % (Table p]) 



The evolutionary status of this star is undetermined. Ryde et al 



( |1999[ ) speculated that IRAS 15194-5115 might be a massive 
(5-8 Mq) star in the last stages of evolution where its low 
'^C/'^'C-ratio is the result of hot bottom burning (HBB). The 
increase of '"^N from the CNO cycle is a signature of HBB 
( Marigo| 200 1 ; Ventura et al. 2002 ), but given the uncertainty in 
the data presented here, this suggestion cannot be confirmed. 



The seven high mass-loss rate carbon stars presented here ex- 
hibit rich spectra at millimetre wavelengths with many molec- 
ular species readily detected. A total of 47 emission lines from 
24 molecular species were detected for the sample stars. The 
mass-loss rate and physical structure of the circumstellar en- 
velope, such as the density and temperature profiles, was care- 
fully estimated based upon a detailed radiative transfer anal- 
ysis of CO. The determination of the mass-loss rate enables 
abundances for the remaining molecular species to be calcu- 
lated. The derived abundances typically agree within a factor 
of five indicating that circumstellar envelopes around carbon 
stars have similar molecular compositions. 

The most striking difference between the abundances are 
reflecting the spread in the '^C/'^'C-ratio of about an order 
of magnitude between the sample stars. Also, the high abun- 
dance of SiO in the envelopes indicates that a shock has passed 
through the gas in the inner parts of the envelope. This is further 
corroborated by the relatively low amounts of CS and possibly 
HCN. 
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The abundances of species that are produced in the outer 
parts of the wind can be reasonably well explained by current 
photochemical models. 
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Appendix A: Spectra 



Fig. A.l. High-resolution spectra of IRAS07454-7 1 12, obtained with the SEST. 



Fig. A.2. ffigh-resolution spectra of IRAS07454-71 12, obtained with the SEST. 



Fig. A.3. High-resolution spectra of CIT 6, obtained at OSO. 



Fig. A.4. Low-resolution specti^a of CIT 6, obtained at OSO. 



Fig. A.5. Low-resolution spectixim of CIT 6, obtained at OSO. 



Fig. A.6. High-resolution spectra of IRAS15082-it808, obtained with the SEST. 



Fig. A.7. High-resolution spectra of IRAS 15082-4808, obtained with the SEST. 



Fig. A.8. Low-resolution spectra of IRAS 15082-4808, obtained with the SEST. 



Fig. A.9. Low-resolution specti-a of IRAS15082-4808, obtained with the SEST. 



Fig. A.IO. High-resolution spectra of AFGL3068, obtained at OSO. 



Fig. A.ll. Low-resolution spectra of AFGL3068, obtained at OSO. 



Fig, A.12, High-resolution spectra of IRC+40540, obtained at OSO. 



Fig. A.13. Low-resolution spectra of IRC +40540, obtained at OSO. 



